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guiding microstructured fibers. While the majority of the field energy in the 
honeycomb photonic bandgap fibers is found to reside in the silica, a sub- 
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may show such fibers particularly interesting for sensor applications, especially 
those based on nonlinear effects or interaction with other structures (e.g. Bragg 
gratings) in the glass. 
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1. Introduction 

Photonic crystal fibers (PCFs), which guide light in a single-material structure by 
coherent scattering from an array of /im-sized airholes (for recent reviews we refer 
to Refs. 111121 and references therein), have in recent years emerged as an attractive 
alternative to conventional optical fibers within the area of nonlinear fiber devices.!^ 
The advantages of the PCFs are firstly that very small mode areas can be obtained 
due to the large refractive index contrast between silica and air, leading to high 
nonlinearity coefficients. Secondly, the PCFs allow for a more flexible tailoring of 
the dispersion properties, which are crucial for many applications. PCFs with zero- 
dispersion wavelengths ranging from 565 to 1550 nm and high nonlinearity coefficients 
have been demonstrated.!^ 

The highly nonlinear PCFs fabricated today are of the index-guiding type,^ in 
which a missing hole in a triangular lattice of airholes defines a high-index core, which 
guides light by total internal reflection. Fibers guiding light in large hollow cores by 
means of the photonic band gap (PEG) effect have also been demonstrated,!^ with the 
intention of obtaining very low losses and nonlinearities. However, an alternative PEG 
fiber design in which a honeycomb airhole lattice is modified by addition of an airhole 
to form a low-index core regiorP has until now not been investigated thoroughly 
although its practical feasibility was demonstrated as early as 1998.!^ The design is 
shown in Fig. d and can be characterised by three parameters: The physical distance 
between nearest-neighbor airholes (commonly denoted the pitch, or A), the diameter 
of the cladding holes, dd, and the diameter of the hole forming the core defect, dc- 
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These fibers guide the majority of the light in silica, just as the conventional index- 
guiding PCFs, and may, therefore, constitute an alternative way of fabricating highly 
nonlinear fibers. The purpose of the present work is to investigate the design depicted 
in Fig. 2] with respect to field energy distribution and nonlinear coefficients. It will 
be shown that nonlinearities comparable to those obtained in index-guiding PCFs 
can be achieved in the honeycomb PBG fibers, while at the same time, a substantial 
fraction of the field energy may be pushed into the airholes. This may make these 
fibers particularly interesting for applications as sensing devices .1^^^^ 

The rest of the paper is organized as follows: In Sect. 2, we describe our theoretical 
methods, and derive a generalized formula for the fiber nonlinearity coefficient, in 
terms of an effective area, which is valid for all field distributions, including those 
where a substantial part of the field energy resides in air. In Sect. 3, we present and 
discuss our numerical results for some selected honeycomb designs, and compare them 
to results for index-guiding PCFs. Finally, Sect. 4 summarizes our conclusions. 

2. Theoretical approach 

We consider a structure, which is uniform along the z-axis while structured in the 
x-y plane. The magnetic field vector, H, may then be written: 



H(r) = e*('^*-^")H(a;,?/), 



(1) 



for a monochromatic wave. The fundamental equation for H is: 



(H, 0H) 
(H,H) ' 



(2) 



0H = V X 




(3) 
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where e(r) = eoer{r) is the dielectric function. For a fixed propagation constant, P, this 
equation may be solved for uj, which can then be regarded as a function of (3 and e{r). 
From the magnetic field vector, the corresponding electric field is straightforwardly 
obtained using Amperes law (SI units are used throughout the paper): 



V X H = iueE (4) 

The effective area is a concept originating in the theory of third-order nonlin- 
earity effects in optical waveguides.^ In a homogeneous material, such as amorphous 
silica, the third-order part of the nonlinear susceptibility gives rise to an amplitude 
dependent shift in the material refractive index: 



n = riQ + An = rii -|- n2 | E p (5) 

where ni is the refractive index of silica in the limit of zero field, and n2 is a nonlinear 
coefficient related to the nonlinear-susceptibility tensor, x^^\ throughP^ 

= ^RexEL- (6) 

In an optical fiber, the change in material index leads to a corresponding change 
in the effective index (riefr = 5)- first-order perturbation theory, the mode-field 
distribution can be considered unchanged by the index perturbation, and the change 
in UJ for a (Kerr-induced) Ae for fixed f3 is from Eq. Q found to be: 



^£oc^(E,Ae,E) 
= 2(H,H) • 
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Usually, the experimental situation is that light is launched at a fixed frequency 
cuq and the Kerr-induced change in refractive index effects a change in the propagation 
constant (3 from /5o, say, to Pq + A/5. A/5 is then determined by: 



lu{(3q + A(3) + Auj = Luo, Lu{(3o) = cuq. (8) 

Here, uj{/3) is the relation between uj and P in the absence of the Kerr effect, and we 
have neglected the change in Acu arising from the shift in (3, assuming that both Acu 
and A/3 are small. We can then obtain A/5 by linear expansion: 



duJ ^ ^ * A ^ Aa; 

—A/3 = -Au =^ A/3 = 9 



a/3 ^ v\ 
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Here, is the waveguide group velocity in the absence of material dispersion effects 
(since we consider propagation at a fixed frequency, the frequency dependence of e 
should not be taken into account when evaluating |^). The change in arising from 
the Kerr effect is then: 



cAoj , . 

^^eff = -nrr- (10) 



For Asr we have: 



Asr = 2^rn2 I E |2 +0 (ti^|E|^) , (11) 
and thereby, neglecting the small n| | E |^ term and using ^/e^ — rii. 



Sqc^ J nin2 \ E l"^ dA Sqc J nin2 \ E \^ dA 
^"^ = Mm) = <(E,D) (12) 
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where, in the last equahty, the inner product of H with itself has been rewritten 
using Eq. (j2|) and (ji}. The integral is taken over the xy-plane. Note that, for a 
microstructured fiber, ni,n2 are position-dependent quantities. 

In the case of ordinary silica fibers, the parameters ni,n2 have little variation 
over the fiber cross section, and it is common practice to express the Kerr-induced 
change in the guided-mode effective index as: 

p 

T? 

A^eff = P-^ (13) 

where P is the power launched into the fiber, = and Aqq is the effective 

mode area. In a microstructured fiber, both rii and n2 can have a strong position 
dependence, and a generalization of Eq. (fT^ is not straightforward. However, in the 
case where ni,n2 are piecewise constant functions, taking on different values over 
the fiber cross section, we can modify the above definition to: 

N P 

Anefr = PE^ (14) 

i=l ^eflf 

where denote the value of in the different sections of the fiber. It can be 
showrP^ that P and Vg are connected by: 

P = J {ExH)JA = t;J(E, D) (15) 
From Eqs. (|T2| . |T5|l we now obtain: 



Ancff = Pl^n^, -2 = P}^n2i " 

i=i (t^°(E,D)) ^=l (E,D,)2 



(16) 



In the last step, we have introduced = Er'E and the effective group index of 
the guided mode, = Note that the integral over | E |^ in each term is restricted 
to the regions with rii = n\. Comparing Eqs. (fT^ and (fTH|l . it can be seen that: 

(nln0)2 I E |4 dA 
_(n\Y (E,D,)2 ^^^^ 



\n'J J,\E-Br\'dA 
In the present paper, we shall only be concerned with the case of pure silica/air 
fibers, in which ni,n2 are equal to 1 and 0, respectively, in the air regions while having 
the values appropriate for silica in the rest of the transverse plane. In this case, the 
nonlinear coefficient will be entirely determined by the effective area relating to the 
silica regions, that is: 

_(n,y (E.D,)' 

The values for ni, ^2 are now understood to be those of pure silica. 

At this point, two comments are in order: Firstly, note that of the two ni factors 
in the denominator of Eq. ()18p. one comes from the definition of in terms of n2, 
whereas the other one comes from the fundamental wave equation. Therefore, if one 
uses a table value of derived with a rii different from that at which experiments 
are done, it may be necessary to use two different ni values in the product. Secondly, 
Eq. ()18|1 differs somewhat from the formula commonly used for Aes , namely: 



(/I 


E |2 dAf 


/ 


1 E 


|4 dA 



= " ' ^ ' (19) 



However, if we assume that all the field energy resides in the silica regions of the 
fiber, and that n° ^ ni we obtain: 



(nX (E,D,)^ _ Unl\^\^dAf ^ J/ | E p 

Thus, the commonly used formula (jl9|) appears as a limiting case of the more 
general result (fTHj) . The approximations leading from Eq. (fTHj) to ()19|) are reasonably 
well fulfilled in standard fibers and in most index-guiding PCFs. This is, however, not 
the case for the fiber designs examined in the present work. In Fig. El we show the 
relative difference between the A^^q definitions in Eqs. (fTH | .(fT ^ defined as: 

AAeff = Aeff - ieff (21) 

with y4efr , Aeff givcu by Eqs. (fT8|) . (fT9|l respectively. The differences have been calcu- 
lated for some fibers designed to have a substantial part of the field energy in air, as 
discussed in the next section. It is evident that differences of 10-20% between the two 
definitions can easily occur. Also, the difference varies in both sign and magnitude 
over the transmission window of the fibers, so that no simple scaling rule between the 
two definitions can be extracted. 

In this work, we solve Eq. ((2)) by expanding the magnetic field and the dielectric 
function in plane waves, using a freely available software package.'^ Eq. Q can be 
used to derive the electric field vector from the magnetic, and the effective area can 
then be found from Eq. (jl8j) . Since it is important for many applications of nonlinear 
effects to be close to a wavelength at which the group velocity dispersion is zero, 
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we have also investigated the dispersion properties of the fibers. In order to scan a 
range of physical pitches efficiently we have used a recently developed perturbative 
approach to the inclusion of material dispersion effects.^ We have found this scheme 
to be both efficient and accurate in the case of silica/air microstructured fibers. 

3. Numerical results and discussion 

In the present investigation, we will primarily focus on structures with lying in 
the interval between 0.3 and 0.8 for both core and cladding holes. Such structures 
are by now routinely fabricated by fusing and drawing hand-stacked capillary tubes 
and rods of silica. However, in order to investigate the limitations of the hole-defect 
honeycomb PEG fibers, we have also studied a design with ^=0.95 and ^ in the 
range between 0.1 and 0.3. Results for f=0.1, 0.2 and 0.3 are shown in Fig. 01 The 
lowest value of A^s is 0.76A^ and is obtained for ^=0.2. The occurence of a minimum 
Aes as a function of core hole diameter can be understood as follows: Since nonlinear 
effects only occur in the silica part of the fiber, reduction of core hole size increase 
the region of integration in the denominator of Eq. (fTH|) . thereby acting to decrease 
the effective area. On the other hand, as is evident from Fig. El reduction of the core 
hole size also decreases the values of ^, where the guided mode becomes localized. 
In other words, the fiber dimensions becomes larger relative to the wavelength of the 
guided mode, and this acts to increase the effective area, when measured relative to 

In the case of index-guiding PCFs, effective areas as low as 1.7 fim^ at a wave- 
length of 1.55 yum have been reported experimentally, corresponding to A^s ~ 0.75A^. 
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The minimal mode area that can be obtained in a sihca rod in air has been proposed 
as a theoretical lower limit on the effective area in silica-based index-guiding fibers, 
and has been found to be 1.48 /xm^ at 1.55 /im (or 0.62A^).'^ Thus, in spite of the 
fact that the guided mode is localized around an airhole in the honeycomb structures 
considered here they are able to obtain mode areas, which are only slightly larger 
than what is possible in index-guiding PCFs. 

In Fig. 13 we show the effective function of wavelength for ^=0.56, 

0.68 and 0.80, and ^ around 0.3-0.5. It can be seen that A^s ~ 1.5 — 3A^ is readily 
obtained, and that the minimal area decreases with increasing size of the cladding 
holes and increases for increasing size of the core hole. 

In many experiments involving nonlinear effects it is important to work at a 
wavelength at which the chromatic dispersion of the fiber is close to zero. In these 
situations, it is the minimal effective area obtainable at a given wavelength, under 
the condition that the dispersion coefficient be zero, which is of interest. Both index- 
and PBG-guiding PCFs can have complicated dispersion curves with several zero- 
dispersion points. Some examples of what can be achieved with the honeycomb design 
for ^=0.68 are shown in Fig. El Following the three curves with dc/dci=0A5 it can be 
seen that the zero-dispersion point can have a discontinuous behaviour as a function 
of A. The curve for dc/dci=0.55 is an example of a fiber with several dispersion zeros. 
To investigate the effective area at the zero-dispersion point we have chosen to focus 
on the longest zero- dispersion wavelength for a given design, since this is where one 
will usually have the smallest effective area relative to the wavelength of the light. 
For the PEG fiber we have, therefore, investigated the location of the longest zero- 
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dispersion wavelength, Aq over a range of physical fiber dimensions from a/3A=1.5 to 
5 /im. In Fig. iniwe report Aq versus the physical pitch, A, and in Fig. [3 the effective 
area at Aq is plotted versus Aq. Broken curves indicate a discontinuity in Aq as a 
function of A. In some cases, the zero- dispersion wavelength can sweep over the same 
frequency several times as A is varied, which is why some of the curves in Fig. [3 are 
multi-valued. Generally, the honeycomb fibers tend to have dispersion zeros falling in 
the range between 0.8 and 2 /im. 

For comparison we have also considered index- guiding PCFs with a triangular 
array of airholes constituting the cladding and a core-defect formed by a missing air 
hole. For air hole diameters d < d* (with d* ~ 0.45A) this class of PCFs is endlessly 
single mod^and by scaling the pitch both large- mode area PCFs as well as small-core 
non-linear PCFs can be formed. In addition to the mode-size the dispersion properties 
may also be engineered. In the top panel of Fig. |Slwe show the zero-dispersion wave- 
length Ao as a function of the pitch A for 5 hole diameters. Depending on the pitch the 
PCF may have none, one, or two (or even three if also considering the near-infrared 
regime) dispersion zeros. For the situation with two dispersion zeros it is seen that the 
second dispersion zero (counting from the short wavelength limit) depends strongly 
on both the pitch and the hole diameter whereas for the first dispersion zero the 
dependence on hole diameter dominates over the much weaker dependence on pitch. 
One of the exciting properties of these PCFs is that by increasing the hole diameter 
the lowest dispersion zero can be moved toward the visible regime. In the lower panel 
of Fig. IHl we show the effective area versus zero-dispersion wavelength. In general we 
find that when the hole size is increased the mode becomes more tightly confined 
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with a smaller effective area. The plot also illustrates a highly desired property; when 
shifting the first dispersion-zero toward the visible the effective area also decreases 
so that the intensity thresholds for various non-linear phenomena also decreases. We 
note that the first dispersion zero may be moved to the visible and the effective area 
may be decreased by further increasing the air hole diameter (so that d > d*), but 
then care must be taken that the PCF remain single-mode near the dispersion zero. 

Comparing Figs. [7| and |H1 it can be seen that the honeycomb PBG fibers inves- 
tigated here do not offer substantially improved flexibility in the tailoring of mode 
areas and dispersion properties compared to the index-guiding PCFs. The PBG de- 
signs with large cladding holes do seem to offer smaller mode areas at the longer values 
of Ao but this could also be obtained in index-guiding PCFs by increasing the size 
of the airholes, thereby going further out of the endlessly single-mode regime. Deter- 
mining the single-mode regime for PBG fibers with large cladding airholes is difficult 
due to the appearance of multiple bandgaps even at relatively long values of |-. For 
the fibers with ^=0.8 we have found that a guided second-order mode is present 
at Ao over most of the A range investigated. For the fibers with smaller values of ^ 
the second-order modes mostly appears at wavelengths shorter than Aq for A-values 
smaller than ~2 /im meaning that a useful range of Aq values without second-order 
mode guidance is available. We have not, however, checked for the presence of guided 
modes in higher-order gaps. 

One interesting feature of the honeycomb PBG design compared to the index- 
guiding PCFs is that a relatively large fraction of the field energy resides in the airholes 
of the fiber. This is illustrated in Fig. IHl where the fraction of electric field energy 
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present in air has been plotted for the fiber designs discussed above. It is evident that 
energy fractions of 10-15% in air are readily obtained, even for holes of moderate size, 
and that the fraction increases with increasing size of the core hole defect. To estimate 
the range of energy distributions accessible, we have investigated some designs, in 
which the core hole defect has been further enlarged. In Fig. EH we show results for 
effective areas and the energy fraction in air for fibers with ^ around 0.6 and varying 
size of the cladding holes. It can be seen that the energy fraction in air can be as large 
as 30%, while still having a range of possible zero-dispersion wavelengths and fairly 
small mode areas. Further increase of d^. does not push appreciably more field energy 
into the airholes, however, the transmission windows and accessible zero-dispersion 
ranges are quickly diminished. The somewhat counterintuitive fact that the that the 
minimum of the effective area curves (i.e., the maximum of the nonlinear coefficient at 
a particular wavelength) occurs approximately at the same values as the maximum 
of the energy fraction in the airholes is due to the fact that the energy fraction in air 
is high at long wavelengths, where the fiber size relative to wavelength is small. 

In Fig. ^2 we show some radial mode profiles obtained by integrating the electric 
field energy density over the angular coordinate in a polar coordinate system around 
the core center. The curves are calculated for fibers with ^ =0.8 and varying size 
of the core hole, and have been normalized so that their radial integral is unity. The 
wavelength has been chosen so as to maximize the fraction of field energy in air. 
Only 5-10% of the field energy is present in the central hole defining the core. This 
shows that the increase in the energy fraction present in air for increasing ^ is not 
so much due to the field energy being pushed into the central core hole, but rather 
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to an increase of the field energy present in the cladding region as is also evident 
from Fig. ^2 The presence of a substantial part of the field energy in air not only 
influences the integrals in Eq. (fT^ . but also makes the group velocity of the guided 
mode deviate substantially from the material refractive index, thus influencing the 
prefactor (^^^ ■ Fig. El this fraction is depicted for the designs with ^=0.8. Its 
importance for a correct evaluation of the effective area in these structures is evident. 

The power fraction in air (which is not completely equivalent to the quantity 
calculated here) for index-guiding PCFs was investigated by Monro and co-worker^ 
who found that large airhole diameters and ^-values were needed to obtain appre- 
ciable power fractions in air. In Fig. ^1 we report the fraction of field energy in air 
for some index-guiding PCFs with the same cladding structure as those discussed 
above but somewhat larger airholes. Comparing Figs. IHl and ^1 with Fig. it is 
evident that the A/ A values needed to push a given fraction of field energy into the 
air region are considerably larger in the index-guiding fibers, calling for smaller values 
of the physical pitch for operation at a given wavelength. This may make the honey- 
comb fibers particularly relevant for evanescent field devices,!^^ such as gas sensors, 
based on interactions with the glass, through, e.g., nonlinear effects or inscribed Bragg 
gratings. 

4. Conclusion 

In conclusion, we have investigated the field energy distribution and nonlinear coef- 
ficients of honeycomb photonic bandgap fibers and compared them to index-guiding 
photonic crystal fibers with a cladding structure consisting of a triangular array of 
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airholes. A generalized concept of effective mode area, which is adequate for the treat- 
ment of fibers with a substantial part of the field energy present in the airholes, has 
been derived for this purpose. While the honeycomb fibers do not seem to offer in- 
creased fiexibility in the design of dispersion properties and mode areas they do offer 
the same possibilities as the index-guiding fibers, at wavelengths above ~1 /xm. In 
addition, the honeycomb fibers have a larger fraction of the field energy present in the 
airholes which may make these fibers particularly interesting for sensor applications 
based on interactions with the glass, through e.g. nonlinear effects and/or inscribed 
Bragg gratings. 
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List of Figure Captions 

Fig. ^ The generic PCF structure investigated in the present work. The core and 
innermost cladding holes are shown along with the defining parameters c/c, dd, and 
A. 

Fig. 121 Relative difference between the effective area definition proposed here, Eq. 
(IT^ . and the commonly used definition, Eq. (fTIH) . 

Fig. 01 Effective area, relative to wavelength, calculated from Eq. ()18|) for a fiber 
with ^=0.95 and various values of the core hole diameter, dc- 

Fig. m Effective area relative to the wavelength of the guided mode for honeycomb 
PBG fibers with ^=0.56 (a), 0.68 (b) and 0.8 (c) and various values of 

Fig. Plots of the chromatic dispersion coefficient, D, in units of ps/nm/km for 
various values of the pitch in a honeycomb PBG fiber with ^=0.68 and ^=0.45. 

Fig. ini Longest zero-dispersion wavelength, Aq, as a function of the physical pitch, 
A. Structures and labeling as in Fig. HI 

Fig. [71 Effective area at the zero- dispersion wavelength, Aq as a function of Aq. 
Structures and labeling as in Fig. HI 

Fig.lHlEffective area and dispersion zeros for index-guiding PCFs. Top panel: Zero- 
dispersion wavelength versus pitch. Lower panel: Effective area at the zero- dispersion 
wavelength as a function of the zero- dispersion wavelength. 

Fig. IHl The fraction of the electric field energy of the guided mode present in the 
airholes. Structures and labeling as in Fig. ^ 

Fig. ^1 Effective area relative to wavelength (a), effective area at the zero- 
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dispersion wavelength (b) and energy fraction in air (c) for some fiber designs in 
which a substantial part of the field energy resides in the airholes. 

Fig. ^2 Radial profile of the electric field energy density, obtained by integration 
over the angular coordinate in a coordinate system with origin at the core center. The 
curves are normalized to have unit radial integrals. The thin vertical lines indicate 
the position of the first ring of cladding airholes. 

Fig. 1121 Variation of the prefactor ^ in Eq. ()18|) with wavelength for a fiber with 
^=0.8 and various values of the core hole diameter, (ic. 

Fig. El Energy fraction in air of index guiding fibers with a cladding structure as 
shown in Fig. |H1 and various airhole diameters d. 
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Fig. 2. LiEGSGAARD 
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Fig. 4. LiEGSGAARD 
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Fig. 6. LiEGSGAARD 
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Fig. 8. LiEGSGAARD 
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Fig. 10. LiEGSGAARD 



29 




Fig. 11. LiEGSGAARD 
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Fig. 12. LiEGSGAARD 
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Fig. 13. LiEGSGAARD 
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